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We reexamine two essential issues within the Peierls-Nabarro model which are critical in obtaining 
accurate values for the Peierls stress. The first issue is related to the sampling scheme of the misfit 
energy across the glide plane and the second one is the effect of relaxation on the Peierls stress. It is 
shown that most of the current applications of the Peierls-Nabarro model do not treat properly the 
two issues and therefore are not able to predict reliable values for the Peierls stress. We argue that 
the double counting scheme for the misfit energy at both sides of the glide plane is physically more 
reasonable, and it can predict more accurate values for the Peierls stress, especially for dislocations 
with equal spacing between alternating atomic rows. We also show the importance for allowing 
atomic relaxation when a dislocation traverses the Peierls barrier, which in turn lowers the Peierls 
stress for narrow dislocations by an order of magnitude. 

(February 7, 2008) 
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The Peierls stress a p is the minimum applied exter- 
nal stress required to move a stationary dislocation ir- 
reversibly, without the assistance from thermal or quan- 
tum lattice vibrations. This fundamental quantity of a 
lattice dislocation was first estimated by Peierls M and 
Nabarro Q using essentially a continuum model, which 
is now referred to as the Peierls-Nabarro (P-N) model. 
The original P-N model has served more as a conceptual 
tool for a qualitative understanding of dislocation core 
properties, rather than providing quantitative estimates 
of these properties, due to the unrealistic sinusoidal force 
law adopted in the model. However, there has been re- 
newed interest in the P-N model recently to calculate 
the Peierls energy and stress thanks to the advent 
of highly accurate ab initio total-energy methods based 
on the density functional theory (DFT). To date, the 
P-N model serves as a link between atomistic and con- 
tinuum approaches, by providing a means to incorporate 
information obtained from the atomistic (ab initio or em- 
pirical) calculations directly into the continuum model. 
This resultant approach can then be applied to problems 
that neither atomistic nor conventional continuum mod- 
els could handle separately. 

After the pioneering work of Peierls and Nabarro, sev- 
eral attempts to improve the original model have been 
made. For example, the three-dimensional dislocation 
profile has replaced the one-dimensional constrained path 
approximation in the original model. This is par- 

ticularly important for treating the dissociation of a per- 
fect dislocation into partials. A fully anisotropic treat- 
ment of the elastic interaction has been implemented into 



the model |T(J. The above extensions have allowed the 
study of straight dislocations with arbitrary orientation 
in an arbitrary glide plane and a medium of arbitrary 
anisotropy. More recently, an averaging procedure has 
been proposed to account for the nonlocal nature of the 
misfit energy [ pd] | . By discretizing the elastic energy, Bu- 
latov and Kaxiras [0 introduced recently a semidiscrete 
variational model which takes into account the impor- 
tant degrees of freedom which participate actively in the 
translation of the dislocation over the Peierls barrier. 

In spite of the apparent importance and success of 
the P-N framework, a number of shortcomings still ex- 
ist in the model. The classic P-N model completely ne- 
glects the relaxation of atoms when the dislocation cen- 
ter moves across the Peierls energy barrier, e.g., the dis- 
location displacement profile stays constant during the 
translation of the dislocation. Thus two approximate 
treatments for the total energy are involved in the clas- 
sic model. Specifically, the variation of the elastic en- 
ergy with respect to the translation vanishes due to the 
fact that the elastic energy is the functional of the over- 
all displacement field of the lattice which are not al- 
lowed to change during the translation of the dislocation, 
hence the elastic energy has no contribution to the Peierls 
stress. On the other hand although the variation of the 
misfit energy with respect to the dislocation translation 
exists due to the discrete lattice nature, no relaxation of 
atoms is allowed to compute the correct (relaxed) lat- 
tice potential resulting from the interaction between the 
moving dislocation and the underlying lattice. In other 
words, in the classic P-N model the dislocation translates 
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rigidly through the lattice. How the two approximations 
affect the description for the dislocation core properties in 
still an open question, therefore it is one of the purposes 
in this work to examine the importance of relaxation in 
determining the Peierls stress and other dislocation core 
properties for both narrow and intermediate dislocations, 
such as those in Al. 

The second important issue addressed in this paper is 
more controversial and subtle. It is related to the misfit 
energy sampling for different crystal lattices, i.e. facing 
lattice (F-type), in which the atoms just above and be- 
low the glide plane face each other, such as the simple 
cubic lattice or alternating lattices (A-type) in which the 
atomic sites on the two sides of the glide plane alternat- 
ingly distribute as in a closed packing lattice, such as 
the fee lattice. In their original work, Peierls, Nabarro 
and others 0J|,|l^] summed the misfit energy indepen- 
dently over the top and bottom half-crystals for both F 
and A-type lattices. We will refer to this scheme as the 
double counting (DC) scheme. However, the DC scheme 
yields the variations of the misfit energy and the lattice 
friction which have a periodicity of 6/2, in contrast with 
the feature of the dislocation barrier which must in gen- 
eral exhibit the periodicity of b. Most people believe 
the unexpected form of the Peierls energy barrier is at- 
tributable to the DC approach in the classic P-N model 
p3| , p^ | and since then almost all numerical and analytical 
work employ exclusively a single counting (SC) scheme 
, e.g., the scheme in which one sums the misfit en- 
ergy only on one side of the glide plane as a function of 
the disregistry (relative displacement) between pairs of 
atomic rows. The argument for using the SC approach is 
based on the assumption of the nearest-neighbor interac- 
tion and on the fact that this scheme gives the correct pe- 
riodicity of b. However, as Wang pointed out recently, the 
undesired periodicity is actually due to an erroneous rep- 
resentation of the atomic positions in the original model 
rather than the DC scheme itself |^,|l(|. In fact, the DC 
scheme can give the correct periodicity for the both F 
and A-type lattices after the error of atomic positions is 
corrected. On the other hand the accuracy of the nearest- 
neighbor interaction approximation in the SC scheme for 
the A-type lattice is still questionable. Thus the second 
aim of this paper is to examine the soundness of the two 
schemes and compare for the first time the misfit energy 
and the Peierls stress calculated from these schemes. Fi- 
nally we present a simple analytical expression for the 
Peierls stress derived for the A-type lattice based on the 
DC scheme. 

In order to study these subtlies within the P-N model, 
we employ the semidiscrcte variational P-N model, be- 
cause of its accuracy and expedition ■ For the study 
of dislocations in Al (intermediate core width) the 7 
surface is obtained by means of ab initio pseudopoten- 
tial total-energy calculations ||. For the case of ex- 
tremely narrow dislocations, we will resort to a sinusoidal 



force law for calculating the misfit energy and the Peierls 
stress. Within the semidiscrete variational approach the 
equilibrium structure of a dislocation is obtained by min- 
imizing the dislocation energy density functional 

U*u = £ \*» + pfpf) + K sP ^ P f] 

hj 

+ £A*7 3 (^£^^(pN°) (D 

i i,l 

with respect to the dislocation densities or displacement 
vectors. Here pj , p^ and p^ are the edge, vertical 
and screw components of the general interplanar dis- 
placement density at the i-th nodal point, and 73 is 
the three-dimensional misfit potential. The correspond- 
ing applied stress components interacting with the py~ , 
pf^ and pf \ are r (1) = a 2 i, t (2) = a 2 2 and r (3) = a 23 , 
respectively. K e and K s are the edge and screw pre- 
logarithmic energy factors. The dislocation density at 
the i-th nodal point is pi — (8% — — a:*— 1), where 

8i and the displacement vector and the coordinate 

of the z-th nodal point (atomic row), respectively, and 

Xij = %(t>i,i-l4>j,j-l + lpi-l,j-l + lpi,j - i'i.j-l ~ lpj,i-l, 
with (f>i j = Xi — Xj, and ipij = \§\ j In \4>ij |. Because 

the displacement vector 8{Xi) is allowed to change dur- 
ing the process of dislocation translation, the Peierls en- 
ergy barrier can be significantly lowered compared to its 
corresponding value from a rigid translation. The re- 
sponse of the dislocation to an applied stress is achieved 
by the minimization of the energy functional with respect 
to pi at a given value of the applied stress, . The 
Peierls stress is then obtained by evaluating the critical 
value of the applied stress rf at which the dislocation 
energy functional fails to be minimized with respect to 
P i through standard conjugate gradient techniques. This 
approach of calculating the Peierls stress is more accurate 
and physical transparent because it captures the nature 
of the Peierls stress as the stress at which the displace- 
ment field of the dislocation undergoes a discontinuous 
transition. 

In Table I, we compare the values of the misfit en- 
ergy and the Peierls stress for various dislocations in Al, 
calculated from the DC and SC schemes. The ab ini- 
tio determined 7 surface has been implemented into the 
semidiscrete model to obtain accurate values for the ener- 
getics and the Peierls stress. All dislocations considered 
here are perfect fee dislocations in (111) plane but with 
different orientations, e.g., different angles between the 
dislocation lines and the Burgers vectors. For all dislo- 
cations in Al, the atomic rows from the two sides of the 
glide plane are situated alternatingly from each other. 
Furthermore it is important to note that, with the excep- 
tion of the 30° and 90° dislocations, the spacing between 
alternating atomic rows for all the other dislocations in 
an fee lattice is not even. We will refer to the 30° and 
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90° dislocations as evenly-spaced dislocations. Two con- 
clusions can be drawn from the comparison in Table I. 
First the SC scheme always underestimates the misfit en- 
ergy and overestimates the Peierls stress compared to the 
DC scheme. Second the agreement for the Peierls stress 
between the two schemes becomes particularly worse for 
the evenly-spaced (30° and 90°) dislocations. Namely, 
while the SC scheme underestimates the misfit energy by 
about 15%, it overestimates the Peierls stress relative to 
its value from the DC scheme by orders of magnitude 
. Thus the even versus uneven spacing between atomic 
rows turns out to be an important issue regarding the 
Peierls stress, and has not been addressed in previous 
studies. 

In order to establish which scheme is superior over the 
other and to understand the different results received 
from two sampling schemes, next we take a closer look at 
the two schemes for the A-type lattice (the two schemes 
are essentially the same for the F-type lattice). In the 
SC scheme, the misfit energy is sampled as a function of 
the relative displacement between pairs of atomic planes 
across the glide plane. Therefore, the SC model consid- 
ers only the nearest-neighbor interaction across the glide 
plane, e.g., the local bonding distortion between the pair 
of atomic rows. Although this model seems to be appli- 
cable to covalently-bonded systems, where the bonding 
across the glide plane is highly localized, it fails to de- 
scribe metallic systems, which have more delocalized elec- 
tronic states. Along the same line of thinking, one can 
understand why the SC model works particularly worse 
for the evenly-spaced (30° and 90°) dislocations in which 
the atomic rows below (above) the glide plane are located 
right in the middle of the atomic rows above (below) the 
glide plane. In this case it is most ambiguous to de- 
fine a local atomic pair, and the second nearest-neighbor 
interaction is as important as the first nearest-neighbor 
interaction. The neglect of higher nearest-neighbors in- 
teractions in the SC scheme gives rise a lower misfit en- 
ergy and a much higher Peierls stress. On the other hand 
the DC scheme treats the misfit energy from the top and 
bottom crystal as separate entities, and the misfit energy 
is defined as a function of the displacement between the 
original position of an atomic row and its final position 
after the introduction of a dislocation. Consequently the 
misfit energy can be summed up independently over the 
top and bottom half-crystals. Two advantages immedi- 
ately emerge from the DC model: (1) Higher-order inter- 
actions are included naturally. In fact the misfit energy is 
exact, by including the contributions from all atomic rows 
across the glide plane. In other words the misfit energy is 
the result from the overall charge density redistribution 
due to the displacements of all atoms on both sides of 
the glide plane, rather than the local bonding distortion 
associated with the SC scheme. (2) The fact that the DC 
approach samples the misfit energy over twice as many 
atomic rows as in the SC method, it reduces the error 



due to the local displacement gradient approximation by 
a factor of two JT^JT^ ] . This improvement is particularly 
important for describing narrow core dislocations where 
the displacement gradients are relatively larger. 

Next we present for the first time a simple and rigorous 
analytical formula of the Peierls stress for dislocations in 
A-type lattice based on the DC scheme. Although this 
formula is derived using the sinusoidal approximation for 
the restoring force, it will provide insight on the effect 
of atomic spacing on the Peierls stress and will allow a 
qualitative understanding of the results in Table I. Fol- 
lowing the SC treatment of Joos and Duesbery [Q, the 
total misfit energy within the DC scheme can be written 
as the sum of misfit energy contributions from the two 
half-crystals 



W{u)= 



ai + a 2 



{j[S(n(ai + a 2 ) - u)\ 



+7[<5(n(ai + a 2 ) + ai - u)]}, 



(2) 



where a\ and a 2 are the spacings for the atomic rows al- 
ternating across the glide plane, which can be different in 
general. As alluded earlier, 6 is the displacement field of 
atomic rows relative their original positions, 7 is the mis- 
fit potential, and u is the dislocation translation distance. 
Introducing a = a± + a 2 and ui = — (ai — u) =u — ai, we 
have 
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W(u,ut)= ^2 -^{^[S(na - u)} + -f[8(na - u/)}} 



n— — OO 



W 1 (u) + W 2 (uf). 



(3) 



Assuming a sinusoidal force law for the restoring force 
F(S(x)) — T max sm(2w6(x))/b, gives a misfit energy func- 
tional r y(S(x)) of the form, 



>y[6(x)] = 



2~ 



-(1 



2wS(x) , 
cos^— ). 



(4) 



The solution of the P-N integro-differential equation gives 
8{x) = b/n arctan(x/C) + b/2, where b is the Burgers vec- 
tor of the dislocation and £ is the half-width. Substi- 
tuting the expressions for 7(2;) and S(x) in Eqn. (3) and 
using the Poisson's summation formula, the misfit energy 
reduces to 



W(yi,y 2 ) = 



Kb 2 sinh 27rr 



1 



8ir cosh 27rr — cos 27ryi 

1 



cosh 27rr — cos 2iry 2 



), 



(5) 



where y\ = u/a, y% = uf/a = {u — a\)/a = y\ — a±/a, and 
r = (/a. if is a constant depending on the elastic prop- 
erties of the dislocations. It should be pointed out that 
the Eqn. (5) is valid only for non-evenly spaced dislo- 
cations due to the requirement imposed by the Poisson's 
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summation formula. It is straightforward however, to 
calculate the Peierls stress for evenly-spaced dislocations 
using a similar approach. To get the Peierls stress one 
needs to calculate the stress associated with the misfit 
energy variation. 



<r(yi,V2) 



1 dW 
b du 



Aosinh27rr 
sin 2-Ky 2 



sin 27r?/i 



(cosh27rr — cos27T7/2) 



(cosh27rr — cos27ryi) 2 

>]• (6) 



Letting t = a%/a and t^h, 



v(yi,t) 



Kb sinh27rr, 



sin 2nyi 



+ 



4a (cosh 27iT — cos 2iry\) 2 
sin27r(yi — t) . 



(cosh27rr - cos27r(yi - i)) 2 ' 



(7) 



For all dislocations studied in Al, we find r = - > 1 and 
this should hold also for not very narrow dislocations. In 
this limit, Eqn. (7) reduces to 



<r{yi,t) = - 



Kb cos irt 
2a cosh 27rr 



sin(27rj/i — Trt). 



(8) 



For a given t, the maximum value of c(yi, t) corresponds 
to the Peierls stress 



Kh 

a p = — e- 27r ^°|cos(7rt)| 
a 



(9) 



For comparison we present below also the formula derived 
by Joos and Duesbery using the SC model [Eqn. (24) in 
Ref. 4] 



Kb 
al 



(10) 



Here al = a . Two results are evident from the sim- 
ple analytical expression derived above. First, Eqn. (8) 
shows that the DC scheme gives the correct periodic- 
ity a. Second, the value of the Peierls stress calculated 
from the DC scheme is always smaller than that from 
the SC scheme. Both these analytical results are consis- 
tent with our numerical results listed in Table 1. For the 
non-evenly spaced dislocations in the fee lattice, t=2/3, 
thus we find that the SC scheme predicts a Peierls stress 
which is twice as large as that obtained from the DC 
scheme using the sinusoidal force law. For the evenly- 
spaced dislocations, it can be easily shown that 



2Kb 



(11) 



and the ratio m(C) of the Peierls stress from the SC (<J p ) 
to the DC (cTp C ) scheme is 

m(C) 



^ C 1 



,2«;/a 



DC 



(12) 



This equation shows that for the evenly spaced dislo- 
cations (30° and 90°), the SC scheme predicts a much 
higher Peierls stress than the DC scheme. It also demon- 
strates that m(Q) increases exponentially with £, yielding 
a larger m(() ratio for the 90° dislocation compared to 
that for the 30°, due to the larger half- width of the 90° 
dislocation Although this simple formula gives the 
correct qualitative trend for a p , it is rather limited in pre- 
dicting reliable Peierls stress for real materials because of 
the sinusoidal force law employed in its derivation. For 
example one problem associated with the sinusoidal force 
law is that it yields no dissociation into partials, which 
can result in several orders of magnitude difference in 
the <7 P if the dissociation is allowed. It needs to be em- 
phasized that in order to test the reliability of different 
models, one needs to compare to direct atomistic calcu- 
lations. The semidiscrete model with the DC scheme has 
been found to be able to predict reliable values for the 
a p for both covalent and metallic systems, e.g., in good 
agreement with the direct atomistic simulation results 
using the same interatomic potentials H||. 

We next examine the effect of atomic relaxation on 
the Peierls stress for different types of dislocations in Al. 
In Table II, we list the a p for different dislocations us- 
ing three different approaches. All the three methods 
are based on the DC scheme. The first (SD+DFT) and 
second (SD+sin) methods employ the semidiscrete model 
but with different restoring forces. Namely the SD+DFT 
method uses a restoring force obtained from density func- 
tional calculations, whereas the second method (SD+sin) 
uses a sinusoidal force law, but with the T max in Eqn. (4) 
adjusted to give the same value for the dislocation half- 
width (C) as in the SD+DFT calculations. The third 
method uses the analytical expressions (Eqns. (9) and 
(11)) for the <j v using the same C as for the first two meth- 
ods. Therefore the comparison between the (SD+DFT) 
and (SD+sin) methods provides information for the ef- 
fect of sinusoidal approximation on a p (both schemes in- 
clude relaxation) , while the comparison between the sec- 
ond and third methods illustrates the effect of atomic 
relaxation. It is interesting to note that the relaxation 
effect is small for all dislocations in Al, whereas the ap- 
proximation of the sinusoidal restoring force is not good, 
as we addressed earlier. 

Recently Schoeck claimed that the relaxation effect on 
the Peierls stress can be canceled out by the opposite 
contributions from the misfit energy and the elastic en- 
ergy pl[ . Thus we need to clarify the cause of the neg- 
ligible relaxation effect in Al, i.e., whether it is due to 
the small Peierls energy barrier or the cancellation ef- 
fect. We have carried out two sets of calculations based 
on the DC scheme and a sinusoidal force law, but r max is 
chosen large enough (1.2 ev/A 3 ) to produce narrow dis- 
locations (£ less than 0.5 A ) in the Al lattice. The first 
set of calculations employs the semidiscrete model which 
allows relaxation while the second set is based on the an- 
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alytical formula o~ p 



3\/3^2 



T max a lKb derived for narrow 



dislocations |J] (no relaxation is included). Note for very 
narrow dislocation the misfit energy is localized within 
one lattice spacing and hence the details of the misfit 
energy sampling (DC vs SC) is not relevant. Listed in 
Table III are the results of the a p for different disloca- 
tions from these two sets of calculations. One can see 
that the relaxation effect is significant (more than one 
order of magnitude) and no canceling effect is observed. 

In summary we have examined two essential issues 
within the Peierls-Nabarro model which are critical to 
obtain accurate values for the Peierls stress. The first 
issue concerns the sampling scheme for the misfit energy. 
We have shown that the DC scheme is more appropriate 
for determining the Peierls stress, especially for evenly- 
spaced dislocations, although most work in the litera- 
ture has done otherwise. An analytical formula is derived 
to help understanding the results received from the nu- 
merical calculations. The second issue addressed in this 
work is the effect of relaxation for various dislocations. 
We have demonstrated the importance of relaxation in 
obtaining reliable Peierls stress for narrow dislocations. 
Comparison of the Peierls stress using ab initio and the 
sinusoidal restoring force illustrates the failure of the si- 
nusoidal approximation. 

The work at the California State University Northridge 
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TABLE I. Comparison of the misfit energy (ev/A) and 
Peierls stress (mev/A ) employing the DC and SC schemes. 
The semidiscrete model with the ab initio 7 surface is used 
for the calculations. 



DC SC 



0° 


misfit energy 


0.0994 


0.0908 




Peierls stress 


1.60 


3.90 


25.3° 


misfit energy 


0.1229 


0.1173 




Peierls stress 


0.04 


0.14 


oU 


misnt energy 


n 1 001 


n 1 no a 
U.1U24 




Peierls stress 


0.33 


3.40 


44.7° 


misfit energy 


0.1370 


0.1280 




Peierls stress 


0.37 


0.84 


60° 


misfit energy 


0.1521 


0.1356 




Peierls stress 


0.61 


1.30 


85.3° 


misfit energy 


0.1706 


0.1630 




Peierls stress 


0.02 


0.07 


90° 


misfit energy 


0.1688 


0.1470 




Peierls stress 


0.02 


2.20 


TABLE II. 


The Peierls stress (mev/A 3 ) 


for dislocations 


in Al determined through three different approaches. The 


DC scheme is 


used in the semidiscrete model (SD) and the 


analytical formula calculations. For the same dislocation all 


the three methods use the same core width in order to give a 


fair comparison. 








SD+DFT SD+sin 


Formula 


0° 


1.60 


0.65 


0.99 


30° 


0.33 


0.02 


0.04 


44.7° 


0.37 


0.04 


0.06 


55.3° 


0.53 


0.07 


0.12 


60° 


0.61 


0.09 


0.14 


66.6° 


0.44 


0.03 


0.05 


TABLE III. 


The Peierls stress (ev 


/A) for dislocations in Al 


lattice calculated from the semidiscrete model (SD) and the 


analytical formula derived only for 


narrow c 


dislocations (see 


text). The same sinusoidal restoring 


force with the magnitude 


7~max — 1.2 CV 


/A 3 is used for both methods. 


The core widths 


for all the dislocations studied here 


are less than 0.5 A. 




SD+sin 




Formula 


0° 


1.2 




18.0 


30° 


0.9 




18.4 


60° 


1.1 




12.9 


90° 


0.7 




13.7 



G 



